UTCCP-P-20 



Four-dimensional Simulation of the Hot Electroweak Phase 
Transition with the SU(2) Gauge-Higgs Model 



Yasumichi Aoki 

Center for Computational Physics, University of Tsukuba, Ibaraki 305, Japan 
(December 1996; revised May 1997) 



Abstract 



We study the finite-temperature electroweak phase transition of the min- 
imal standard model within the four-dimensional SU(2) gauge-Higgs model. 
Monte Carlo simulations are performed for intermediate values of the Higgs 
boson mass in the range 50 ^ Mh ^ lOOGeV on a lattice with the temporal 
size Nt = 2. The order of the transition is systematically examined using 
finite-size scaling methods. Behavior of the interface tension and the latent 
heat for an increasing Higgs boson mass is also investigated. Our results 
suggest that the first-order transition terminates around Mh ~ 80GeV. 
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I. INTRODUCTION 



The possibility that the baryon number asymmetry of the universe was generated in the 
course of the electroweak phase transition in the early universe has been much discussed. 
One of the conditions necessary for a realization of this possibility is the presence of out-of- 
equilibrium process during the transition. If the transition is of first order with a sufficient 
strength, super cooling could occur and the universe could have been driven out of equilib- 
rium. This raises the question whether the electroweak phase transition undergoes a strong 
enough first-order transition for a realistic Higgs boson mass. 

In the minimal standard model the essential features of the finite-temperature elec- 
troweak phase transition are expected to be controlled by the SU(2) gauge-Higgs sector. 
Hence finite-temperature studies of the SU(2) gauge-Higgs model represent a first step in 
understanding the electroweak transition. In the past two approaches have been taken for a 
non-perturbative analysis of the transition in this model by means of Monte Carlo simula- 
tions (for reviews, see Refs. [|1],3). One approach is to employ the original four- dimensional 
model. The other is to treat a three-dimensional model derived by a dimensional reduction 
of the original model in time. In the three-dimensional approach many studies have been 
made covering a wide range of the Higgs boson mass 35 < Mh < ISOGeV 0-01 • In par- 
ticular, it has recently been reported that the phase transition is absent for Mh ^ 80GeV 
On the other hand, studies with the four- dimensional model have been less 



extensive, and systematic analyses have been limited to relatively light Higgs boson masses 
of Mh ^ 20 — 50GeV |[T0|-[T2|. Since the Higgs boson mass is experimentally bounded by 



Mh > 64GeV [T^ , it is important to extend the study of the four-dimensional model toward 
heavier Higgs boson masses. 

In this article, we report results of our simulation of the four-dimensional model for the 
Higgs boson mass in the range 50 ^ Mh ^ lOOGeV. Our aim is a systematic analysis of 
the strength of the transition, which is known to be of strong first order for light Higgs 
boson, as the Higgs boson mass Mh increases. For this purpose finite-size scaling analyses 
of susceptibility and the Binder cumulant are carried out. We also study latent heat and 
interface tension in order to examine the strength of the transition in terms of physical 
quantities. 

In Sec. |T| we describe some details of our simulation. Monte Carlo time histories of 
observables are discussed in Sec. iHIl Results of finite-size scaling analyses are presented 



in Sec. [V. In Sec. and VI we report measurement of latent heat and interface tension. 



Our results are summarized in Sec. |VI1| . In Appendix we describe our procedure to estimate 
the Higgs boson mass Mh and the critical temperature for our simulation points from 
published results. 

II. SIMULATION 

We employ the standard lattice action given by 

^ = E ( E f Trf/„,,. + Y: «:Tr(<|.t f/„,,<|.„+^) -pi- X{pl - 1)4 , (1) 
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where the radial part p„ is defined by the decomposition of the complex 2x2 matrix Higgs 
field $„ = Pn<^n, Pn > 0; ttn G SU(2). All of our simulations are made for the temporal 
extent Nt = 2. We set the gauge coupling j3 = l/g"^ = 8. Simulations are made for 6 
values of the scalar self-coupling A, choosing the hopping parameter k in the vicinity of the 
transition point. The parameter values of all of our runs are listed in Table |. Also listed 
in the table are estimates of the zero-temperature Higgs boson mass Mh and the critical 
temperature Tc at the transition point of Nt = 2 lattice. These estimates are obtained by 
interpolating available data for Mh and at A = 0.0001, 0.0003, 0.0005(Mh = 47GeV) 
|n|,|g, 0.0017235(Mj:^ = 85GeV) and 0.023705 §. The procedure of the interpolation is 
described in Appendix. 

For each value of A runs are made on an A^^ x 2 lattice with A^^ = 8, 12, 16, 24, 32, and in 
addition with A^, = 40 for A = 0.001, 0.0017235, 0.003(Mh = 66,85,108GeV ). Gauge and 



scalar fields are updated with a combination of the heat bath and overrelaxation fTS 



algorithms in the ratio displayed in Table which has been reported to have the fastest 
decorrelation time in terms of number of sweeps in Ref. We make 10^ iterations of 

the combined updates for each parameter point except for the run with A = 0.0005, A^, = 
32, K = 0.128865, where we make 2 x 10^ iterations. 

Observables are measured at each iteration. Our analyses are mainly made in terms of 
the quantities defined by 



P 
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mNt ^ 3 



-Tr($tt/„,$ 



A. = 



n "J 1=1 



pi- 



rn 

(3) 
(4) 



where Lg is the spatial part of the scalar hopping term in the action, A^ represents the 
angular part of Lg and the squared Higgs condensate. 



III. MONTE CARLO TIME HISTORY 

A qualitative test for a first-order phase transition is provided by Monte Carlo time 
histories of observables and their histograms. In Fig. |l| we show the time history of and 
its histogram for the largest volume for each M^- For the lightest Higgs boson Mh = 47GeV 
(Fig. |l|(a)), a clear fiip-fiop behavior is seen, which is refiected in the double peak structure 
of the histogram of Lg. If we increase the Higgs boson mass, the fiip-fiop behavior becomes 
milder, indicating a weakening of the first-order transition. It can be seen, however, up to 
Mh = 66GeV (Fig. |l](d)). A double peak structure of the histogram also persists up to this 
value of Mh- These observations provide qualitative indication for a first-order transition 
up to Mh = 66GeV. 

On the other hand, for the heavier Higgs boson of Mh = 85 and 108GeV, we cannot 
discern a fiip-fiop behavior nor a double peak structure in the histogram. We have analyzed 
the histogram varying the hopping parameter k around the value of the simulation through 
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the reweighting technique [jT6[. For Mh = 85GeV (Fig. |T|(e)), we found that the histogram 



changes to a narrow flat peak at k = 0.129966, while only a sharp single peak structure has 
been seen for = lOSGeV (Fig. |l|(f)). This suggests that a first-order transition is either 
absent, or extremely weak at these Higgs boson masses. 



IV. FINITE-SIZE SCALING ANALYSIS 

Let us consider the susceptibility of defined by 

XA. = V ((A.^) - (A.)2) , (5) 

with V = Ng the spatial volume. On a finite lattice this quantity has a peak at the pseudo 
critical point Kc{V). If the transition is of first order, the maximum value of the peak should 
increase linearly in volume, and the pseudo critical point should converge to the value at 
infinite volume linearly in 1/V [1^. On the other hand, the maximum value should remain 
constant if there is no phase transition. 

We calculate the maximum value and the position of the peak with the standard reweight- 
ing technique |T^. Errors are estimated by a jackknife analysis using 1000-20000 sweeps as 



the bin size, which is chosen from stability of the magnitude of error as a function of bin 
size and varies with coupling parameters and lattice sizes. 

In Fig. results for the maximum value of the susceptibility of A^ is shown as a function of 
the volume V. For 47 < Mh < 66GeV the increase of the maximum value is consistent with 
a linear behavior, quantitatively supporting a first-order transition for this region of Higgs 
boson mass. In contrast, we observe a very flat volume dependence for 85 < Mh < lOSGeV, 
albeit the maximum value is increasing slowly in the range of volume used here. 

Another useful indicator of a first-order transition is the Binder cumulant defined by 

R -1 ^^^^'^ rfil 

^'' = ^- 3 (A?)^- 

For the case of a first-order phase transition the value of deviates from 2/3 in the infinite 
volume limit, while it converges to 2/3 otherwise. 

In Fig. ^ we plot the valley depth of the Binder cumulant of A^ as a function of the inverse 
volume 1/V obtained by a reweighting procedure similar to that for the susceptibility xa^- 
Lines are linear fits to the largest three volumes for each Mh- For 47 < Mh < 66GeV the 
value extrapolated to the infinite volume -B^^ clearly deviates from 2/3, providing additional 
evidence for a first-order transition. For 85 < Mh < lOSGeV the deviation decreases by 
an order of magnitude, although still finite within the error. In Fig. ^ the deviation of B'^^ 
from 2/3 is plotted as a function of the Higgs boson mass. It can be seen that the strength 
of the first-order transition rapidly diminishes toward larger Mh- 

Our results for the susceptibility and the Binder cumulant clearly show that the transition 
is of first order for 47 < Mh < 66GeV. It is also clear that the transition, if first order, 
is a very weak one at Mh = 85GeV and lOSGeV. It is possible that the transition turns 
into a crossover for this range of Mh- Data for larger volumes are needed, however, for a 
conclusive analysis on this point. 
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In Fig. 1^ we summarize our results for the pseudo critical point k,c{V) obtained from the 
susceptibility and the Binder cumulant. In Table |T| we summarize estimates of the critical 
hopping parameter in the infinite volume obtained by a linear extrapolation of k,c{V) in 1/V 
using results of largest three volumes for each A. 



V. LATENT HEAT 



We calculate the latent heat with the equation |]T8 



Ae ~ -M^KA(p2). (7) 

where A means the difference of (p^) in the pure symmetric phase and in the pure broken 
phase. It has been shown P|JT^ that the latent heat calculated with this equation is in good 
agreement with that obtained with the energy operator. 

In order to calculate the expectation value (p^) in the pure phases, we separate the 
ensemble of configurations generated in a run into two sub-ensembles of pure phases. The 
division is made by inspecting the Monte Carlo time history of A^. In order to avoid 
contaminations from transition stages, iterations in those stages are removed following the 
method employed in Ref. [|l^ . 



A source of systematic error in the procedure above arises from the fact that the point of 
K at which the run is performed does not exactly coincide with the pseudo critical point Kc 
which we estimate from the susceptibility xAs- order to treat this problem we make two 
independent runs, both in the region of metastability, such that Kc is sandwiched between 
the hopping parameter of the two runs. Then we interpolate results of the two runs for (p^) 
with a linear function of n to Kc- The error of the latent heat is calculated from those of 
(p^) and Kc- 

The result of our analysis is shown in Fig. M. At the lightest Higgs boson mass of 



Mff = 47GeV our result is in good agreement with that of Ref. []TT| obtained with the 
energy operator. The latent heat decrease apparently linearly with an increasing Higgs 
boson mass Mh- A linear extrapolation of our data is consistent with a vanishing of the 
latent heat at Mh ~ 80GeV. 



VI. INTERFACE TENSION 



The interface tension a provides another indicator of the strength of the first-order 
transition. We calculate this quantity with the Binder's histogram method ||20|. Let Pmax 
and Pmin be the peak and valley height of the distribution of A^, reweighted such that the 
two peaks have an equal height. Define ay = —{N^/2N'^) \ia{Pmin/Pmax)- For spatially cubic 
lattices used in our simulations, finite-size formula for the true interface tension a is given 

by m 



c - - In A^, 
4 



(8) 



where c is a constant independent of A^^. Making a two parameter fit of ay obtained for 
the largest three volumes for Mh = 47, 52, 57GeV or two volumes for Mh = 66GeV we find 
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(t/Tj? shown in Fig. |^. The interface tension rapidly decrease with an increasing Higgs boson 
mass and seems to vanish around Mh ~ 80GeV as in the case of the latent heat. 



VII. SUMMARY 

Our finite-size scaling study establishes a first-order transition for 47 < Mh < 66GeV. 
For larger Higgs boson masses a rapid weakening of the transition makes it difficult to 
draw a definitive conclusion on the order within the range of lattice volumes employed in 
our simulation. However, combining finite-size data with results for the latent heat and 
the interface tension, our four- dimensional study suggests that the first-order transition 
terminates around Mh ~ 80GeV in the = 2 SU(2) gauge-Higgs model. This is consistent 
with the results of recent finite-size scaling studies carried out in the dimensionally reduced 
three-dimensional model @J^. 
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APPENDIX A: 

In Table | we listed the zero temperature Higgs boson mass Mh and the critical tem- 
perature Tc for each A at our critical hopping parameter Kc on the Nt = 2 lattice. Here we 
describe how we estimated those values. 



In Table IV we summarize published results for Mh and at five different values of A re- 



ported in Refs. H JTT|Jl^ . Two values of A employed in our simulation (A = 0.0005, 0.0017235 



) are among those in Table |IV|. For these cases we adopt the published results for Mh and 
Tc. For the other four values of A we estimate Mh and Tc by an interpolation of results in 
Table 0. 

At the tree revel the mass ratio Mh/Mw is proportional to \/~X/k (the gauge coupling 
is fixed to /5 = 8 in our simulations). Plotting five values of Mh in Table |rV| against 
\/X/k we find that they exhibit a smooth increase suitable for an interpolation. A cubic 
polynomial of VX/k yields a good fit with x^/dof = 0.10. We then estimate the values of 
MHatX = 0.000625, 0.00075, 0.001, 0.003 from the value of the fitted curve at TA/k, where 
Kc is the critical point determined from the susceptibility of A^ listed in Table |TTl[ One can 
also make a fit with a quadratic polynomial (x^/dof = 0.05). The difference of values of 
Mh from the two fits may be used as a measure of the systematic error of our estimation 
procedure. However, the difference is quite small ( 0.03 — 0.3%) compared with the statistical 
error (2 — 6%), and we adopt the statistical error in our estimate. 

The results for the critical temperature Tc has a sharp bend for large -\/A/ k. Hence this 
variable is not suitable for interpolation of Tc. We find, however, that Tc as a function of 
Mh is smoothly increasing, and a cubic polynomial fit goes well (x^/dof = 0.5). We use 
the fitted curve and the values of Mh estimated above to find the values of T^. Fit with a 



6 



quadratic polynomial (x^/dof = 0.9) is used to estimate the systematic error, which turns 
out to be comparable to the statistical one. Those errors are combined into the final value 
of error shown in Table |I[ 
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TABLES 



TABLE I. Run parameters of simulation. For each parameter 10^ sweeps are made except for 
the run marked by f (2 x 10^ sweeps). Data generated are used for analysis of susceptibility (x), 
Binder cumulant(i?), latent heat(Ae) and interface tension(cj). For our procedure to estimate the 
Higgs boson mass Mh and the critical temperature Tc, see Appendix. 



A Mg(GeV) TcjGeV) k use 

0.0005 47(2) 94(1) 8 0.128950 x 

12 0.128900 x,B 

16 0.128883 X,B,a 

24 0.128866 X,B,a 

32 0.128860 B 

32 0.128862 Ae 

t 32 0.128865 X,Ae,o- 

0.000625 52(1) 100(3) 8 0.129110 x 

12 0.129036 x,B 

16 0.129004 X,B,a 

24 0.128986 X,B,a 

32 0.128983 5,Ae 

32 0.128987 X,Ae,(7 

0.00075 57(1) 107(3) 8 0.129243 x 

12 0.129158 x,B 

16 0.129126 X,B,a 

24 0.129103 X,B,a 

32 0.129098 B 

32 0.129102 X,Ae,(7 

32 0.129106 Ae 
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TABLE! {Continued) 





Te(GeV) 


Ns 




use 


0.001 66(2) 


119(3) 


8 


0.129476 


Y 
A 






12 


0.129407 


X,B 






16 


0.129349 


X,B 






24 


0.129330 


X,B 






32 


0.129328 


X,B,a 






40 


0.129327 


X,B,Ae,a 






40 


0.129330 


Ae 


0.0017235 85(3) 


139(2) 


8 


0.130200 


X 






12 


0.130080 


X,B 






16 


0.130026 


X,B 






24 


0.129980 


X,B 






32 


0.129966 


X,B 






40 


0.129968 


X,B 


003 108('6'l 


157(10) 


8 


131350 


v 

A 






12 


0.131050 


B 






12 


0.131200 


X 






16 


0.131111 


X,B 






24 


0.131065 


X,B 






32 


0.131054 


X,B 






40 


0.131042 


X,B 



TABLE II. Number of iteration of each updating step in one sweep. 



u 


heat bath 


overrelaxation 

U a p 


1 


1 


3 3 1 


TABLE III. Critical hopping parameter Kc for each A. 


A 


0.0005 0.000625 0.00075 0.001 0.0017235 0.003 


l^c {BaJ 


0.1288607(23) 0.1289812(11) 0.1290981(14) 0.1293262(11) 0.1299633(29) 0.1310360(40) 
0.1288584(13) 0.1289822(25) 0.1290974(14) 0.1293266(11) 0.1299622(27) 0.1310318(34) 
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TABLE IV. Summary of zero temperature Higgs boson mass Mh and critical temperature Tc 
at the transition point of A'^t = 2 lattice. 



A K MniGeY) Tc(GeV) reference 



0.0001 0.1283 18(1) 38(1) Q 

0.0003 0.12865 35(1) 71(1) 

0.0005 0.12885 47(2) 94(1) 

0.0017235 0.13 85(3) 139(2) 

0.023705 0.145 196(13) 185(4) 



11 



FIGURES 




(a) ?i=0.0005 



N3=32 



j" 1.5 



K=0. 128865 



III 



sweep (x1 ) 



2.5 
2.0 
1.5 
1.0 
0.5 



(b) X=0.000625 
N =32 



K=0. 128987 




10 8 6 4 2 





10 8 6 4 2 



sweep (x1 ) 



2.5 
2.0 
f 1.5 
1.0 I 
0.5 



(c) ^=0.00075 



N,=32 



K=0.129102 



sweep (xlO ) 



10 8 6 4 




FIG. 1. Monte Carlo time history of the spatial hopping term Lg and its histogram, 
at every 50th sweep are displayed for history. Histograms are made with full ensemble. 
Mh = 47GeV, (b) Mh = 52GeV and (c) Mh = 57GeV. 
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j" 1.5 
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(d) ^=0.001 



N3=40 



K=0. 129327 




sweep (xlO ) 



2.5 
2.0 
f 1.5 
1.0 
0.5 



sweep (x1 ) 



2.5 
2.0 
f 1.5 
1.0 
0.5 




10 8 6 4 2 



(e) ^=0.001 7235 




: N3=40 








K=0. 129968 





10 8 6 4 2 



(f) ^=0.003 




N,=40 








K=0.131042 






10 8 6 4 2 



sweep (xlO ) 

FIG. 1. (Continued) (d) Mh = 66GeV, (e) Mh = 85GeV and (f) Mh = lOSGeV. 
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FIG. 2. Maximum height of the susceptibihty of as a function of the volume. Lines are 
guides for eyes. 




FIG. 3. Valley depth of Binder cumulant of as a function of the inverse volume. Filled 
symbols show the infinite volume limit. 
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FIG. 4. Valley depth of at infinite volume as a function of Higgs boson 
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FIG. 5. Pseudo critical hopping parameter Kc defined by the maximum point of XAs or the 
minimum point of . 
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FIG. 6. Latent heat as a function of Higgs boson mass. Open symbol is from Ref. and is 
shifted shghtly left for visualization. 
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FIG. 7. Interface tension as a function of Higgs boson mass. Legends are the same as in Fig. 
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